University of Manchester a Fortran 77 Code for Plotting Stability Regions Department of Mathematics a Fortran 77 Code for Plotting Stability Regions
نویسندگان
چکیده
The standard technique for obtaining the stability regions of numerical methods for ordinary diierential equations (ODEs) and delay diierential equations (DDEs) is the boundary-locus algorithm (BLA). However, in the case of the DDE y 0 (t) = y(t) + y(t ?) for & 2 R, the BLA often fails to map out the stability region correctly. In this paper we give a FORTRAN 77 listing of an alternative stability boundary algorithm, which can be used for both implicit and explicit numerical methods applied to both ODEs and DDEs.
منابع مشابه
A Fortran Code for Null Geodesic Solutions in the Lemâıtre-Tolman-Bondi Spacetime
This paper describes the Fortran 77 code SIMU, version 1.1, designed for numerical simulations of observational relations along the past null geodesic in the LemâıtreTolman-Bondi (LTB) spacetime. SIMU aims at finding scale invariant solutions of the average density, but due to its full modularity it can be easily adapted to any application which requires LTB’s null geodesic solutions. In versio...
متن کاملA New Two Dimensional Model for Pollutant Transport in Ajichai River
Accurate prediction of pollution control and environmental protection need a good understanding of pollutant dynamics. Numerical model techniques are important apparatus in this research area. So a 2500 line FORTRAN 95 version code was conducted in which using approximate Riemann solver, couples the shallow water and pollution transport agents in two dimensions by the aid of unstructured meshes...
متن کاملFortran 77 Subroutines for Computing the Eigenvalues of Hamiltonian Matrices II
This article describes Fortran 77 subroutines for computing eigenvalues and invariant subspaces of Hamiltonian and skew-Hamiltonian matrices. The implemented algorithms are based on orthogonal symplectic decompositions, implying numerical backward stability as well as symmetry preservation for the computed eigenvalues. These algorithms are supplemented with balancing and block algorithms, which...
متن کاملExperiences with Porting the Parallel Cfd Code N3s to Fortran 90
We first present the difficulties we have overcome using the CF90 parallelizer-compiler to obtain a parallel version of the N3S code which is similar to that obtained using the CF77 parallelizer-compiler. We also discuss about the stability of the results produced after this porting. We finally describe the difficulties encountered in controlling the memory consumption when moving from the Fort...
متن کاملOn the convergence of the homotopy analysis method to solve the system of partial differential equations
One of the efficient and powerful schemes to solve linear and nonlinear equations is homotopy analysis method (HAM). In this work, we obtain the approximate solution of a system of partial differential equations (PDEs) by means of HAM. For this purpose, we develop the concept of HAM for a system of PDEs as a matrix form. Then, we prove the convergence theorem and apply the proposed method to fi...
متن کامل